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The wavefunction of a massless fermion consists of two chiralities, left-handed and right-handed, 
which are eigenstates of the chiral operator. The theory of weak interactions of elementally particle 
physics is not symmetric about the two chiralities, and such a symmetry breaking theory is referred 
to as a chiral gauge theory. The chiral gauge theory can be applied to the massless Dirac particles 
of graphene. In this paper we show within the framework of the chiral gauge theory for graphene 
that a topological soliton exists near the boundary of a graphene nanoribbon in the presence of a 
strain. This soliton is a zero-energy state connecting two chiralities and is an elementally excitation 
transporting a pseudospin. The soliton should be observable by means of a scanning tunneling 
microscopy experiment. 



For a massless fermion, the left-handed and right- 
handed chiralities are a good quantum number and the 
two chirality eigenstates evolve independently according 
to the Weyl equations. One chirality state goes into the 
other chirality state under a change in parity. The weak 
interactions in elementary particle physics act differently 
on the left-handed and right-handed states, which results 
in well-known phenomena, such as the parity violation 
for nuclear /3 decayji The weak force is described by a 
gauge field. In general, a gauge field which has a differ- 
ent (the same) sign of the coupling for the left-handed 
and right-handed chiralities is called an axial (a vector) 
gauge field. ■^ In the presence of an axial component, the 
interaction between a gauge field and a fermion can be 
asymmetric for the two chiralities. For example, in the 
case of the weak interactions for neutrinos, only the left- 
handed chirality couples with a gauge field and the theory 
is generally known as a chiral gauge theory. 

A chiral gauge theory framework can be applied to 
graphene. The energy band structure for the electrons 
in graphene^iA has a structure similar to the massless 
fermion, in which the dynamics of the electrons near 
the two Fermi points called the K and K' points in the 
two dimensional Brillouin zone is governed by the Weyl 
equations.^ Because the K and K' points are related 
to each other under parity, two energy states near the 
K and K' points correspond to right-handed and left- 



handed chiralities, respectively. The spin for a fermion 
corresponds to a pseudospin for graphene which is ex- 
pressed by a two-component wavefunction for the A and 
B sublattices of a hexagonal lattice,^ The correspond- 
ing pseudo-magnetic field for the pseudospin is given by 
an axial gauge field that is induced by a deformation of 
the lattice in graphene. ^"^ The electronic properties of 
a graphene are thus described as a chiral gauge theoryi^ 
An important point here is that the axial gauge field in 
graphene has different signs for the coupling constants 
about the two chiralities while the conventional electro- 
magnetic (vector) gauge field does not. 

In a chiral gauge theory, the chiral symmetry break- 
ing and the resultant mixing of chiralities are of prime 
importance. In elementally particle physics, this symme- 
try breaking relates to the origin of the mass of a fermion 
and the experimental investigations into the mass of neu- 
trinos are in progress. Since graphene is described by a 
chiral gauge theory, a chirality mixing phenomenon in 
graphene is a matter of interest. In this paper, we show 
that a graphene nanoribbon which is a graphene with a 
finite width having two edges at the both sides fiS^i- has 
a chirality mixed soliton solution when applying strain to 
a graphene nanoribbon. Two symmetric edge structures, 
that is, armchair and zigzag edges are shown in Fig.[T] It 
is known that the spatially localized electronic states, the 
edge states, appear near the zigzag edget^— A chirality 
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FIG. 1: Structures of a polyacetylene and a graphene edge. 
a, Two possible isomers trans- and cis-polyacetylene. b, Two 
principal edge structures: zigzag and armchair edges. H de- 
notes a hydrogen atom, and carbon atoms are divided into A 
(•) and B (o) atoms. 



mixed soliton consists of two edge states belonging to dif- 
ferent chiralities, and it is a natural extension of the con- 
cept of the topological soliton in trans-polyacetylene.— "— 

(Definition of gauge fields) 

First we review the chiral gauge theory of graphene;^ A 
lattice deformation in graphene gives rise to a change of 
the nearest-neighbor hopping integral from the average 
value, —7, as —■^ + 6ja{'r), where a (= 1, 2, 3) denotes the 
direction of a bond as shown in Fig. [2^. We define the 
axial gauge A(r) = {Ax{r),Ay{r)) by Sja{r) 



as^-^ 
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vpA^ir) = (571 (r) - - {(572 (r) -I- (573(r)} , 
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{S-f2{r) - S-ysir)} : 



where vp is the Fermi velocity. The direction of the vector 
A(r) is perpendicular to that of the C-C bond with a 
modified hopping integral, as shown in Fig. [2}3. The 
effective Hamiltonian for a deformed graphene is written 
by a 4x4 matrix as^ 



iJ*(r) = Up 



(P- 



A(r) - eA<"" 
CTx0*(r) 



(r)) 



r 



(P-A(r)- 



3A°™(r)) 



*K(r) 
*K'(r) 



(2) 



where the field (j){r) relates to A(r) as vp(j){r) = {Ax{r) + 
iAy{r))e~'^^^^^ in which kp is the Fermi wave vector of 
the K point, and A'""(r) is an electromagnetic gauge 
field. Here, cr = {ax,Oy) [a' — {~ax,ay)] are the Pauli 
matrices which operate on the two-component spinors 
^K (r) and ^I'k' (r) for the pseudospin. We use the units 
Vp = 1 and h = 1, and thus the momentum opera- 
tor becomes p ~ — iV. A lattice deformation does not 
break time-reversal symmetry, which appears as the dif- 
ferent signs in front of the field A(r) for the two chirali- 
ties, while the electromagnetic gauge field A''™(r) breaks 
time-reversal symmetry and has the same sign for the 
K and K' points. A (A°™) is an axial (a vector) gauge 
field. ^ Like the case of A®™(r), the field strength of A(r) 
defined as Bz{r) = dxAy{r) — dyAx{r), plays a fundamen- 
tal role in discussing topological solitons and edge states, 
as we will show below. 

It is straightforward to show using equation ([T]) that 
the field (f) behaves as a position-independent interac- 
tion for the Kekule distortion,^ and then equation ([2]) is 
equivalent to the Dirac equation with a mass in four- 
dimensional space-time without the z-component [pz — 
0] (see Appendix \^ . Though the main concern of this 
paper is a chirality mixing due to a local mass (^(r), let 
us begin by considering the massless limit (j){r) = and 
examining the chirality eigenstate ^k (right-handed chi- 
rality) using the 2x2 Hamiltonian, H{r) = (T-(p + A(r)). 

(Topology of the gauge field) 



In Fig. [2)2 the double bond represents a shrinking of 
the C-C bond and the single bond denotes the absence 
of deformation. The a phase is defined as the bond- 
ing structure for the case of {S^i,Sj2,Sj3) = (0,^7,0), 
while the /3 phase is the case of {S^i , Sj2, ^73) = (0, 0, Sj). 
From equation ([T]), the corresponding A fields for the a 
and [3 phases, A+ and A_, are given, respectively, by 
A± = (— (^7/2, ±-\/3(57/2). For the skeleton of a trans- 
polyacetylene shown between the closed dashed lines of 
Fig. [5):, it is well-known that a topological soliton ap- 
pears when the configuration has a domain wall (a kink) , 
that is, when the f3 phase changes into the a phase at 
some position along the x-axisi^"— The gauge field for 
such a domain wall configuration for a zigzag nanoribbon 
is written as 



M{x) = (cj.,Ay(a;)), 



(3) 



where Cx 



-(57/2, Ayix) 



X <^ — ^, and Ay{x) 



{ay ^ \/3(57/2) when 



when X ^ ^. Here, ^ (> 0) de- 



notes the width of a kink (see Fig. [^b) . In addition, the 
gauge field which describes the edge structure is given 
by A2. This A2 comes from the fact that the C-C 
bonds at the zigzag edge are cut.^^ This cutting is rep- 
resented by (^71, (572, (573) = (7,0,0) at the edge, and 
A2 = (7,0). Since there are two zigzag edges aX y — yu 
and y — yi in the zigzag nanoribbon (without a domain 
wall), A2(j/) = {Ax{y),0) has a value only for y = yu 
and y = yi (the edge location), and otherwise Ax{y) = 0. 
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FIG. 2: Representing a lattice deformation in terms of the 
axial (deformation-induced) gauge field, a, A lattice deforma- 
tion is defined by {Syi,S'^2, ^73). b. The direction of vector A 
is perpendicular to the lattice deformation. The directions of 
the arrows are for the case of a positive S'y. c. The configura- 
tion of the axial gauge field A for a trans zigzag nanoribbon. 
The two distinct bonding structures, a phase (A+) and /3 
phase (A_), are combined together to form a a domain wall 
(a kink). The y-component of the field Ai(a;) changes its sign 
at a: = 0, which represents a kink structure. The zigzag edges 
are represented by the field A2(y). 



The total gauge field for a trans zigzag nanoribbon is 
given by the sum of Ai(x) in equation ([3]) and A2{y) as 
Ai(a;) + A2(y) = {cx + Ax{y),Ay{x)). As a result, the 
(K point) Hamiltonian is given by 

H{r) = a^iPx +Cx + Ax{y)) + ay{py + Ay{x)). (4) 

(Zero-energy solution of H{r)) 

Here we assume that the energy eigenstates of H{r) 
in equation ^ has the form of e*P="^*p^(a;, y)|(T), where 
Px is the quantum number and \a) denotes the spinor 
eigenstate. The energy eigenequation is rewritten as 

{o-xipx + Dx + Ax(y)) +<Jy{py + Ay (x))} ^ p^ {x , y)\a) 



E^p^{x,y)\a) 



(5) 



where D^ = Px + Cx- We decompose this eigenequation 
into two parts by putting ^p^ (a:, y) = ip{x)ip{y) and E = 
El + E2 as 

WxPx + <JyAy{x)} ip{x)\a) = Eiip{x)\a), (6) 

{ax{Dx + Ax{y)) + aypy} ^{y)\a) ^ E2^iy)\a), (7) 

In general, the spinor eigenstate of the first equation can 
not be identical to that of the second one. However, in 
the special case that Ei — E2 — 0, the spinor eigen- 
states of these equations can be the same. It is because 
that -ff(r) commutes with cr^ for the zero-energy state, 
[iJ(r),cr^]_e*P-^*^="(r)|cr) = 0, and that the spinor 



eigenstate can be taken as the eigen spinor of az defined 

as (Tz|ct±) = ±|tT±), where 



l^+) 



Thus, the corresponding zero-energy states are pseu- 
dospin polarized states, namely, the amplitude appears 
only one of the two sublattices. In the following, we show 
that equations (O and ^ give, respectively, the topo- 
logical soliton^^"— i^ and the edge statesi^i^^ From these 
zero-energy states for equations dH) and ([7]), a general 
zero-energy solution for equation ([5]) can be constructed. 

(Topological soliton) 

Let us obtain the zero-energy soliton for equation ^ . 
When El = 0, the eigenequation is represented as 
WxPx + o-yAy{x)}'iJj±{x)\a±) = 0. We have two solu- 
tions. 



(8) 



i^±{x) — iVexp I ± / Ay{x)dx 



where A^ is a normalization constant. When we use a 
trial function Ay{x) = Oj, tanh(a;/^), we get ^'±(2;) = 
Ncos\v^'^''^{x/£,)M Hence, when ay > (kink), only i/j^ 
is selected, while when Uy < (anti-kink), only -0+ is 
selected. The significance of a single zero-energy state is 
that the particle-hole symmetric partner is given by itself, 
which leads to the result that a soliton has no charge but 
has spin l/2i^i22 The sign of Uy corresponds to the sign 
of the field strength as Bz{x) = Uy / {^ cosh^ {x / ^)) . The 
sign of the Bz field is essential to a rule for obtaining 
the normalizable solution. This is easy to understand 
by noting that the square of ^(r) is given by H{r)^ — 
(p -I- A(r))^ -I- Bz{r)(Tz, which gives a positive coupling 
for +Bz{r)az- Because H{yY = for a zero-energy state 
and (p+A(r))^ is always a positive value, the zero-energy 
state needs to satisfy +Bz{r)<Jz < 0, so that a positive 
Bz (> 0) selects |cr_) (or 0_) and a negative Bz (< 0) 
selects \a+) (or ■(/;+). 

(Edge states) 

The derivation of the zero-energy edge states from 
equation (O is given in Ref. [26|. For the case oi Dx < 0, 
there are degenerate zero-energy states given by 



oD^\y-y^\ 



V.-(y)|a_)=e^^|y-y'l|a_), 



(9) 



where |Z?j;|^^ is the localization length. As shown in 
Fig. ^, at the upper edge located at y = yu, the 
Axijj) field increases abruptly when y approaches y„ 
{y < yu)- Therefore, the corresponding i?^ field [Bz{y) — 
~dyAx{y)] is pointing toward the negative z-axis there. 
Hence, only the \a+) state can appear near the upper 
edge. In contrast, at the lower edge {y — yi), the Ax{y) 
field decreases abruptly as y moves away from yi {y > yi). 
Therefore, the corresponding field strength is positive 
there, and only the \cr-) state is selected near the lower 
edge. 



(Soliton-edge state) 

A zero-energy solution of equation ([5]) is constructed 
by the product of the topological soliton V'- (x) of equa- 
tion ^ and the edge state (p-{y)\a-) of equation ^ 
as 



'^pJx,y)W-) =e 



- J" Ay(x)dx^D^\y-yi\ 



(10) 



This new state is localized not only near the lower zigzag 
edge but also near the kink. Since a kink satisfies 
Bz{x) > 0, this state ^^^(a;, j/)|ct_) is the solution to 
equation ([Sj . If there is an anti-kink with Bz (x) < at 
X = 0, another zero-energy state given by 



^lix,y)\a+) = e+r ^«(-)'^-e^==l^-^"l|a+), 



(11) 



is the solution. This state is localized near another zigzag 
edge and is also localized near the anti-kink. In addi- 
tion to these zero-energy solutions of the K point Hamil- 
tonian, there are zero-energy solutions of the K' point 
Hamiltonian. Let the solutions for the K' point be of 
the form of ^_p^(a;,y)|tT) = tp' {x)ip' {y)\<7) . The energy 
eigenequation for the K' point Hamiltonian, cr' • (p — 
A(r)), leads to a pair of energy eigenequations: 

{crxPx + cryAy{x)}'ijy{x)\a) = ~Ei'ijy{x)\a), 
{(^x{Dx + Aj^ijj)) + (TyPy) Lp' {y)\(T) ^ E2(p'{y)\<j). 

These eigenequations are the same as those given in equa- 
tions (O and ([7]) (except for the unimportant sign change 
of ii'i). As a result, the solutions to these equations 
are the same as equations ((TU)) and (jTT]) . We thus have 
two zero-energy solutions originating from the K and K' 
points for a given p^- This number of zero-energy states 
for a ribbon is different from a single zero-energy state 
for a polyacetylene chain. This difference is attributed to 
the fact that the soliton for a polyacetylene chain results 
from a chirality (or an intervalley) mixing. 

(Chirality mixing) 

The zero-energy solutions given by equations (jlOp and 
dill) have been obtained on the assumption that a chiral- 
ity mixing between the K and K' points can be neglected. 
However, translational symmetry along the a;-axis is bro- 
ken due to the presence of a kink (or anti-kink) and a 
kink itself causes a mixing of chiralities. In this case, 
the eigenfunction may be written as a linear combina- 
tion of ^p^ for the K point and ^_p^ for the K' point. 
Their mixing is determined by the mass term which is 
expressed by means of valleyspin Ta {a = 0,1,2,3) as 
a, {riRe[0(r)] - T2lm[0(r)]}. By putting *k = e'^'^^^K 
and ^K' = e+*'''^^K' into equation ([2]) for a zigzag rib- 
bon, we obtain the equations for a general case: 

{t3 {a^xPx + (TyAyix)) - T2e'^'^^'''=a^Ay{x)} ^ = ^i*. 

To Wxicx -k + A^{y)) + aypy} * == £2^^, 

(12) 

where 6k = k^ — k^ The last term on the left-hand 
side of the first equation of equation ((T^ shows that the 



effective domain wall profile for the mixing term is an os- 
cillating function of a;, in contrast to a smooth function 
of the intravalley mixing term for the second term. It 
is straightforward to find a zero-energy solution of equa- 
tion dH]) as ^^ = U{x)il}±{x)tf±{y)\(j±), where U(x) 
is a matrix for valleyspin which satisfies the equation 
d^U{x) - Tie"^'^^''''Ay{x)Uix) = 0. Due to the chirality 
mixing, the actual wave function of a zero-energy state in 
a zigzag nanoribbon can be complicated. One example 
of the wave function is shown in Fig. [3^. 

(Soliton in polyacetylene) 

Equation (|12p can be solved analytically for the spe- 
cial case oi 6k — {k = fcp). In this case, we obtain 
simultaneous differential equations: 



o-x {yPx - T2Ay{x)^ ^ = Ei^, 
T3 { yp^r + <TyAy{x)^ * = E2^, 
To {cTxicx - kp + Ax{y)) + (Typy} * = £3^. 



(13) 



The first term gives rise to a chirality mixing. For 
a zero-energy solution of the first term, the spinor 
eigenstate should be the eigen spinor of ri defined 
as ti|t±) = ±|r±), which shows a strong chirality 
mixing. The zero-energy solutions for the first two 
terms are given by ip±{x)\a±) (^ |Tip), where ip±{x) = 
N exp (± J^ 2Ay{x)dx^ . This state is a valleyspin unpo- 
larized state and also a pseudospin polarized state, and 
these properties are consistent with those of the topo- 
logical soliton in polyacetylene<^ The third equation in 
equation (|13|) describes the edge state having the shortest 
localization length since the localization length is given 
by \cx — kF\~^ which vanishes in the continuum limit. 
The resulting zero-energy solution of equation (fT3)) given 
by v^(r) = '0±(a;)((9±(j/)|(T±)(g)|Tzp) corresponds to Fig. |3)d 
which reproduces a topological soliton in polyacetylene at 
the zigzag edge sites (see Fig. |3h for comparison). Note 
that the soliton can move along the zigzag edge and the 
soliton has a mass. In the case of polyacetylene, the soli- 
ton mass is estimated to be around 6me , where me is the 
mass of the free electron. ^^ In the case of the ribbon, we 
obtain 65{W/£,)me, where W denotes the ribbon width. 
This result reproduces the soliton mass in polyacetylene 
when W = a and ^ = 10a, where a is the lattice constant. 

To further elucidate the effect of the edge on the soli- 
ton, we consider the solitons of an armchair tube, a 
metallic zigzag tube, and a metallic armchair ribbon in 
Appendices |B] and O We show that the chirality mixing 
is negligible for the zero-energy states in these tubes. A 
metallic armchair ribbon produces chirality mixed soli- 
tons when there is a domain wall. The solitons are not 
localized near the edge since there is no edge states near 
the armchair edge. This feature is in contrast to that of 
the soliton in a zigzag nanoribbon. See Appendices [B] 
and [Cj for more details. 

(Discussion) 

We can use equation ([T]) for a lattice deformation in- 
duced by a strain. Let u(r) = {ux{r),yy{r)) be the dis- 
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FIG. 3: Wave function patterns for zero-energy soliton. a, An 
example of the wave function pattern. The solid and empty 
circles represent the phases (+ or — ) of the wave function, and 
the diameter of each circle is proportional to the amplitude. 
We use 57 = 0.27 and the kink profile of tanh(a;/^) with 
^ — 2A. b, An example of the zero-energy soliton in a zigzag 
ribbon. Because the wave function of this example appears 
only at the edge sites, this state is identical to the topological 
soliton in a trans-polyacetylene shown in c. 



shape graphene nanoribbon caused by an acoustic shear 
deformation given by w^; = and Uy{x) = w ln[cosh(x/^)] 
with u = ^(tty/g), can reproduce the gauge field rep- 
resenting a bond alternation (a domain wall) in zigzag 
ribbons. Because g ~ 7)^ u is smaller than 1^ by a factor 
of S^/j. This shows that a domain wall can be realized 
by a strain of ^ 10%.'^^ Similarly, a strain produces a 
localized soliton in a metallic armchair nanoribbon (see 
Appendix \C\\ . On the other hand, the ribbon does not 
support the edge states without the strain. A pseudospin 
polarized wavefunction pattern that is spatially localized 
near the bottom of a "V" shape graphene nanoribbon is 
the indication of a chirality mixed soliton. Note that a 
strain makes it possible to observe the soliton by means 
of a scanning tunneling microscopy (STM) experiment, 
in contrast to that the STM is unable to detect a soliton 
in polyacetylene since the soliton is moving. Moreover, it 
was suggested recently by Guinea et a/.-^ that a uniform 
Bz field may be realized in graphene by a strain-induced 
lattice deformation, which is an interesting consequence. 
If this is the case, it is expected that the Landau level ap- 
pears only for one chirality and the other chirality decou- 
ples from the gauge fields in the presence of a magnetic 
field which eliminates B^ for one chirality. Then the chi- 
ral symmetry in graphene is maximally broken, and this 
situation is similar to the case of weak interactions in 
elementally particle physics. 
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where g is the electron-phonon coupling. An interesting 
consequence of this is that the field configurations which 
are equivalent to A± may be realized when an appropri- 
ate strain is applied to a sample. For example, a "V" 



Appendix A: Original Dirac Hamiltonian 



The original Dirac Hamiltonian is written as 



J 



H^{v) 



(p + ^(r) 
m 



-V(r)) 



-^•(p-^(r)+V(r)) 
I 



*R(r) 
*L(r) 



where m is the mass of fermion. The electronic Hamil- 
tonian for graphene corresponds to the case in which 
^R -^ ^K, *L -^ CTai^K', and TO -^ cJ){y). The vector 
gauge field V(r) and axial gauge field -4(r) correspond to 



eA''™(r) and A(r), respectively. The third component 
such as Pz is assumed to be zero when we identify the 
original Dirac equation (in 3-1-1 dimensional space-time) 
with the effective Hamiltonian for graphene (in 2 + 1 di- 



mensional space-time). 



Appendix B: Solitons in armchair nanotubes 

Here we consider the solitons in armchair nan- 
otubes. The K point Hamihonian is given by remov- 
ing Ax{y) from equation (5). By putting "^p^{x,y) — 
g-iD^K^j-^.'jgiPHy j]2to the energy eigenequation (5), we 
obtain {d^T {Py + Ay{x))} i^±{x) = for the zero- 
energy state. It fohows that the function '4'±ix) contains 
the exponential function ex.p{zLpyx), so that either ip+i^) 
with pj, = or '0_ (x) with py = can be a normalizable 
solution. The momentum Py is quantized by a periodic 
boundary condition around the tube's axis, and a zero- 
momentum state Py — satisfies the boundary condition 
for any armchair nanotube<^ The solution with py — is 
a topological soliton. From equation (2), we obtain the 
chirality mixing term as a^ {TiKe[(f>{x)] — T2lin[4>{x)]}, 
where (j){x) = iAy{x)e~'^^''^^ . This mixing term is small 
because a smooth function Ay (x) of x is multiplied by a 
rapid oscillating function e"^*'^''^. Moreover, the chiral- 
ity mixing term does not cause a first order energy shift 
since the unperturbed states il^±{x)\a±) are pseudospin 
polarized states satisfying {a±\ax\<J±) = 0. For these 
reasons the chirality mixing is negligible in the case of an 
armchair nanotube. 

Note that the zero-energy solitons in an armchair nan- 
otube obtained above are distinct from the topological 
soliton in polyacetylene. The chirality mixing term is ir- 
relevant to the solitons in armchair nanotubes, while it 
is relevant to the topological soliton in polyacetylene. 



Appendix C: Solitons in zigzag nanotubes and 
armchair ribbons 



Let us examine solitons in a zigzag nanotube and an 
armchair ribbon. The existence of a zero-energy topolog- 
ical soliton in a zigzag tube requires two factors: the A 
field topology and the presence of the Dirac singularity. 
The A field topology can be understood by noting that 
the basic unit of structure is a cis-polyacetylene for which 
the two phases shown in Fig. 2^ can be considered. ^^ The 
a phase is defined by {6ji , 6j2, S^a) = (0, ^7, 6j), and the 
(3 phase is {Sji,6j2, ^73) — (Sj, 0, 0). From equation (1), 
the corresponding gauge fields for the a and (3 phases, A+ 
and A_, are given, respectively, by A± = (=F(57,0) [see 
Fig. |4)d]. a domain wall kink is represented by Ai(y) = 
{A^{y),0) with A^{y) = -Sjtanh{y/Q. By assuming 
that the wave function is of the form oie'^Pyyip{x)(p{y)\a), 
we have a pair of the eigenequations from the K point 
Hamiltonian as 

{(T^A^iy) + aypy} ip{y)\a} = Eiip{y)\(T), 
WxPx + crypy}ip{x)\a) = E2ip{x)\a). 

The first equation possesses a zero-energy topological 
soliton. Therefore, when there is a zero-energy state 
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FIG. 4: Soliton in an armchair nanoribbon. a. The wave func- 
tion of a topological soliton in a metallic armchair nanoribbon 
obtained from a tight-binding model, b, The two phases A+ 
and A_ are separated by a domain wall kink distortion rep- 
resented by the shaded region. 



for the second equation, the K point Hamiltonian may 
possess a mixed zero-energy solution. The state with 
Px = and Py ~ 0, i.e., the state at the Dirac singular- 
ity, can satisfy the second equation with E2 = 0. Since 
Px is quantized by a periodic boundary condition around 
the tube's axis, this state with vanishing wave vector 
exists only for "metallic" zigzag tubes. '^'^ For "semicon- 
ducting" zigzag tubes, the quantized Px misses the Dirac 
singularity, and therefore such a zero-energy topologi- 
cal soliton does not exist. Thus, only the presence of 
a non-vanishing B^ field strength does not necessarily 
result in the presence of a zero-energy state. In addi- 
tion to a domain wall, the Dirac singularity is rather 
essential for the presence of a zero-energy state. Note 
that a non-topological excitation, a polaron, may ex- 
ist even in "semiconducting" zigzag tubes near a bound 
kink-antikink pairi^ 

The localization pattern of a topological soliton is sen- 
sitive to the lattice structure of the edge of a nanoribbon. 
To illustrate this, we show the wave function of a topo- 
logical soliton in a "metallic" armchair nanoribbon in 
Fig. 2^. The soliton is extended along the kink, which is 
contrasted with the localized feature of the wave function 
of a zero-energy state in a zigzag nanoribbon shown in 
Fig. 3a (in the text). This difference is a consequence of 
the fact that unrolling a zigzag tube can be represented 
by a strong intervalley mixing term (/)(r) at the armchair 
edge, and that this field (/)(r) does not destroy the Dirac 
singularity."^^ As a result, a topological soliton appears in 
a "metallic" armchair ribbon, as illustrated in Fig.|3h- It 
is interesting to note that unrolling a "metallic" zigzag 
tube does result in a "semiconducting" armchair ribbon. 
This implies that a topological soliton in a "metallic" 
zigzag tube disappears when the tube is unrolled since 
the Dirac singularity also disappears then. 
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